function output=heaviside_fun(X)
t=0;
for i=1:length(X)
   if X(i)>0
    t=t+1;
   end
end
% heaviside(x>=0)=1 &  heaviside(x<0)=0 
if t==length(i)
     output=1;
 else
    output=0;
 end
